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This paper is a continuation of our work on edge-flames in premixed combustion. An edge-flame is a 
two-dimensional structure constructed from a one-dimensional configuration that has two stable solutions 
(bistable equilibrium ). Edge -flames can display wavelike behavior, advancing as ignition fronts or retreating 
as failure waves. Here we consider two one-dimensional configurations: twin deflagrations in a straining 
flow generated by the counterflow of fresh streams of mixture; and a single deflagration subject to radiation 
losses. The edge-flames constructed from the first configuration have positive or negative speeds, according 
to the value of the strain rate. But our numerical solutions strongly suggest that only positive speeds 
(corresponding to ignition fronts) can exist for the second configuration. We show that this phenomenon 
j can also occur in diffusion flames when the Lewis numbers are small. And we discuss the asymptotics of 
die one-dimensional twin deflagration configuration, an overlooked problem from the 70s. 


. V Introduction 

edge-flame can be defined, roughly speaking, 
a flame sheet with an edge. There is a growing 
afcure on edge-flames in no n-p re mixed combus- 
(g.g., Ref. [1]) reflecting their important role. A 
spreading over a fuel bed, solid or liquid, will 
an edge; edge flames are an important part of 
structure of the combustion field that occurs in 
burning of heterogeneous solid propellants [2]; 
y are a noticeable characteristic of candle flames 
T microgravity conditions [3]; and they must ex- 
’wSenever a hole is tom in a flame-sheet by tur- 
ent eddies, so that their behavior is relevant to 
jproblem of lifted turbulent diffusion flames 
Important theoretical work has been pio- 
ly Dold and coworkers [6,7]. 

-flames in premixed combustion have been 
well studied. Indeed, there appears to be only a 
theoretical treatment [8] and a single explicit 
mental study [9]. Old experiments evidence 
their existence [10] is noted in Ref. [8]. The pres- 
is a sequel to Ref. [8] and examines pre- 
edge-flames in two hitherto unexamined con- 
ions. 

^nple framework in which edge-flames can be 
flnicted proceeds in the following fashion: There 
pertain one-dimensional combustion systems, 
dent on a spatial variable x, say, for which 
*re multiple solutions. Of particular interest 
stable systems for which there are three solu- 
te stable and one unstable, the counterflow 


diffusion flame being a well-known example. Sup- 
pose the weakest (strongest) of the stable solutions 
has a temperature distribution T t (x) [T 2 (x)]. 7\(x) 
will often correspond to a quenched state, or a close 
approximation thereof, with values close to some 
background or supply temperature Tf. The maxi- 
mum value of T 2 (x) will be close to a flame-tem- 
perature (e.g., the adiabatic flame-temperature or 
the Burke-Schumann flame-temperature). .An un- 
steady two-dimensional combustion field can then 
be defined by an initial-value problem in the x — z 
plane, where the boundary conditions in x are those 
for the one-dimensional problem, and 


T(x, z, t) — » T 2 (x), T.(x) as z 


4-oo — c 


If T 2 (x) is associated with a single thin reaction zone 
(flame-sheet) and 7\(x) is associated with negligible 
reaction, an edge-flame is defined and the edge- 
flame structure either moves in the direction of de- 
creasing z, in which case we call it an ignition front, 
or moves in the direction of increasing z, in which 
case we call it a failure wave. Here we continue the 
discussion of ignition fronts and failure waves in the 
premixed context. We discuss the situation where 
the one-dimensional problem is defined by a twin- 
flame counterflow, and we also discuss a deflagration 
in which multiple solutions arise because of radiation 
losses. Some related results on diffusion flames are 
briefly discussed, as is the asymptotic description of 
quenching for the twin-flame problem. As the latter 
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Fig. 1. Variations of maximum temperature with Dam- 
kohler number for twin counterflow flames: E, exact; El, 
first-order approximation of exact solution; A, solution of 
the autonomous equation A10; AL (A19), 1 term in q; A 2, 
(A19), 2 terms in q ; B, (A20). 


will only be of interest to the asymptotics commu- 
nity, it is relegated to the Appendix. 

This work and Ref. [8] are parts of the dissertation 
of Vedarajan [11]. 


The Twin-Flame Counterflow Problem 

Consider a symmetric counterflow of fresh mix- 
ture that, in general, supports twin flames. This con- 
figuration has long been studied, both theoretically 
and experimentally. 

A simple model suitable for our purposes starts 
with the equation 

-paxdJdx(C p T, Y) = XCj'd 2 /dx 2 (C p T , Y) 

+ BY e-v™ (Q y - 1) (2) 

where a is the rate of strain and we have assumed 
that the Lewis number is 1. The supply conditions 


1x1 -> oo T Tfi Y -> Y f (3) 

If AT a 3 QYf/C p is used as a reference temperature, 
and 

lc * •jMpCpO. 

is used as a reference length, a single nondimen- 
sional equation can be deduced, namely, 

—xT x = + D( 1 + T f - T)e ~ 0/T , (4) 

D = B/pa, 6 = E/RAT a 

Because of the symmetry, it is sufficient to solve this 
equation in x > 0 with the boundary conditions 

r r ( 0 ) = 0 , T(«) = Tf (5) 

Note that within a context that resolves the cold- 
boundarv difficulty (e.g., a cutoff temperature), the 


system of equations 4 and 5 has a stable quenched 
solution in which T = Tf everywhere. In addition, it 
is well known that for D greater than some minimum 
value D nim , there are two solutions. These can be 
characterized by the maximum value of T (T max ), and 
a representative response, obtained numerically, is 
shown as a solid line in Fig. 1 (curve E). Here Tf = 
0.2, T a = 1 + Tf = 1.2, and 0 = 16, values adopted 
throughout the paper. The upper branch corre- 
sponds to stable solutions, and it is these, along with 
the quenched solutions, that are the key ingredients 
of the unsteady two-dimensional problem descnbed 
in general terms in the Introduction. We are par- 
ticularly concerned with solutions for values of D 
close to D min . 

As an aside, not part of the main thrust of our 
discussion, we note that there have been asvmptotic 
treatments of the system of equations 2 and 3 in lieu 
of a numerical strategy, but there are ingredients 
that have not been discussed before, and we de- 
scribe these in the Appendix. 


Failure Waves and Ignition Waves 

In this section, we use the strategy described in 
the Introduction to define a two-dimensional un- 
steady problem with T>(x) defined by a point on the 
upper branch of Fig. 1 (curve E) and 7\(x) ^ Tf, the 
quenched state. The appropriate generalization of 
equation 4 is 

T t - xT r = T a + 7L + D(T a - T]e-^ (6) 

now accounting for diffusion in the ~ direction as 
well as in the x direction. 

Note that as Id — > °°, where T zz — > 0, equation 6 
reduces to the one-dimensional form (equation 4), 
and the solution T 2 (x) can be assigned as z — > + 00 , 
the solution T x (x) s 7% as z -> - (Note the remark 
following equation 5.) Boundary conditions in x are 
those for the one-dimensional problem, namely, 
equation 5. And initial conditions, whose precise 
form is not important, are defined by a simple 
smooth interpolation between T i and 7Y 

Figure 2a shows an initial temperature profile. 
Note that at z = - 5, it is T\ m Tf , and at z = + 15, 
it is T 2 (x), the one-dimensional twin-flame configu- 
ration. The Damkohler number is 4.9 X 10' and the 
flame sheets are merged together, because we are 
close to D min . Figure 2b shows the profile at a later 
time, and it is clear that the structure is retreating, 
corresponding to a failure wave. 

In view of the double-flame structure, it might be 
argued that the use of the rubric “edge -flame” is not 
appropriate, but the behavior revealed here has its 
counterpart in the edge-flames discussed in Ref. [8], 
where the single-flame counterflow problem is ex- 
amined (fresh versus inert counterflow). 

If a is decreased (D increased), it is possible to 
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FlC. 2. Temperature profiles of a failure wave D = 4.9 X 10": («;) t = 0: ( b ) t - 4.5. 
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get an ignition front, and Fig. 3 shows reaction-rate 
contours in such a case. Steepening of the leading 
(lower edge) portion of the flame is quite apparent. 
Here (and also for the failure wave) the front is of 
unchanging shape following the passage of initial 
transients. 


By performing a number of calculations for dif- 
ferent values of D, it is possible to construct a graph 
showing variations of the front speed V with D (Fig. 
4). Here 

V = UL S' 1 (7) 

where (z)+ is the dimensional front speed and S is 
the adiabatic flame speed defined by the system. The 
quenching value of D (D, nm ) is 4.802 X 10\ and for 
values of D between D mm and D {) (9.776 X 10”), the 
front speed is negative. But for values of D greater 
than D 0 , the front speed is positive. The significance 
of results such as these, both in the present context 
and in the single-flame context of Ref. [8], is that 
flames can experience two distinct types of quench- 
ing: global quenching, arising when D is reduced 
below the value D, nin ; and quenching if part of the 
flame has been destroyed and D inill < D < D 0 so 
that a failure wave enlarges the region of destruction. 

Results such as these and the one -dimensional an- 
alog discussed in Ref. [8] suggest that the basic in- 
gredients we have included here — multivaluedness, 
quenching — will always lead to the dichotomy of 
positive and negative speeds for the two-dimensional 
structure. That this is not the case is shown by the 
example discussed in the next section. 


Multivaluedness Due to Radiation Losses 

Consider the equation for a plane deflagration, ra- 
diation losses included, 

MC p T x = XT „ + BQYe E/RT - q" (8) 

q'" = 10 5 MWP, lhn {Y f - Y)(T - T f ) W/m 3 

where MW is the mean molecular weight, P. ltm is the 
total pressure in atmospheres, and T is the tempera- 
ture in degrees Kelvin. We suppose that Y represents 
methane and the radiating species are C0 2 and H 2 O t 
whose concentrations are proportional to (Yf — Y). 
There is a corresponding equation for Y (cf. equation 
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FlG. 5. Variations of flame temperature with equivalence 
ratio for a deflagration with radiation losses. 

2). The expression for q” is derived bv first con- 
structing an accurate representation of the tempera- 
ture dependence of the band radiation for each spe- 
cies and then fitting this with a straight line, a 
surprisingly reasonable approximation up to a tem- 
perature of 1800 K. All other parameters are chosen 
so that the model provides a reasonable approxi- 
mation of methane/air flames in the neighborhood 
of the lean limit. The goals are only qualitative, but 
we do not wish to be led astray by gross quantitative 
inaccuracies. 

Now the system of which equation 8 is a part de- 
fines, in addition to the quenched solution T = Ty, 
Y - Y f> dual solutions for values of Yy greater than 



Fig. 7. Temperature contours for an ignition front with 
radiation losses, <t> — 0.53556. Values are ( 1.5, 2.0. 2.5, 3.0. 
3.5, 4.0|. 


the radiation-defined inflammability limit Yy =* 
0.02945 (equivalence ratio <f> = 0.5355): See Fig. 5, 
the upper branch of which corresponds to stable so- 
lutions. Thus we can construct a two-dimensional 
unsteady combustion field (an edge-flame), and Fig. 
6 shows an ignition wave for <t> = 0.53556, a value 
close to the limit. Temperature and reaction-rate 
contours are shown in Figs. 7 and 8. The wave is 
quite thick because the temperature decay occurs on 
a scale that is much larger than the preheat- 
zone thickness. Figure 9a shows variations of the 
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Fig. 6. An ignition front in a flame with radiation losses, <f> = 0.53556. 
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08. Reaction-rate contours for ail ignition front with 
losses. <t> = 0.53556. Values are |0.04. 0.1. 0.3. 
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with (f>, results which strongly suggest 
* speeds do not occur. Thus the limit 
ratio is 0.5355, and yet at an equivalence 
.53556, the edge speed is positive with a 
toproximately 0.3555. Of course, it is not 
to reduce <f> arbitrarily close to the limit in 
number of calculations, but it seems unlikelv, 
P* 9a is examined, that anything but positive 
I can be achieved. 


also a characteristic of existing results when the un- 
derlying one-dimensional flame is a diffusion flame 
(e.g.. Ref. [1]). 

On the other hand, if a hole is tom in a flame 
subject to radiation losses, the hole will heal no mat- 
ter how close the equivalence ratio is to the limit 
value. It is natural to wonder if robust flames of this 
nature can occur in the counterflow configuration if 
the Lewis number is different from 1, and although 
we have no results of this nature for premixed 
flames, we have uncovered robust edges in the case 
of diffusion flames. Figure 9b shows variations of 
edge speed versus Damkohler number when the un- 
derlying one-dimensional flame is defined by a sim- 
ple symmetric counterflow configuration (1:1 stoi- 
chiometrv, both Lewis numbers equal to 0.3), and 
negative speeds are not obtained. However, this 
might not be a commonly realizable phenomenon. 
In examining edge-flames constructed from the S- 
shaped response of the Kirkbv-Schmitz configura- 
tion [12] (flux conditions for the fuel applied to one 
boundary, Dirichiet conditions for the oxidizer ap- 
plied at the other), for a variety of parameter choices, 
we have failed to identify robust edges even for 
Lewis numbers as small as 0.2. 

Clearly, the results presented here and in Ref. [8] 
could have relevance to the behavior of turbulent 
premixed flames in the laminar flamelet regime, par- 
ticularly where the response of turbulent flame 
speed to turbulent intensity "bends.” 


jjj^* ' Conclusions 

paper, we have examined two examples of 
lines (but note the remark following equation 
: can arise when the reactants are premixed, 
putts for the symmetric (twin-flame) counter- 
irifieuration are similar to those of the single- 
pcnfiguration examined in Ref. [8]. These can 
iup&rized by the observation that if a hole is 
iTluch a flame, it will get larger or smaller, 
on the value of the straining rate. This is 


Appendix 

Asymptotic Solution of the One- Dimensional 
Counterflow Problem 

Here we briefly discuss the system of equations 4 
and 5 in the asymptotic limit 0 — > there are in- 

gredients that have not been discussed before. In the 
limit, reaction is confined to a flame sheet located at 
x = x*. and if x* ^ 0, we have, to first order. 
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Fig. 9. Non -negative edge speeds: 
(a) deflagration with radiation losses 
and ( b ) symmetric diffusion counter- 
flow configuration, Le\ = Ley = 
0.3. 
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0 < x < £» T = 1 + Tj = T tn } - 0; (Al) 


= ln~ l 5 -2 


x > X* T - Tf = I - Y/Yf 


= erfc(x/^2)/erfc(x^/ ^5) 


defined for S~ 2 < 7C/2. This branch, with infinite 
slope at S' 2 = rtf 2, is drawn in Fig. 1 (curve El) 
along with the upper branch T* - T„ (f* = 0) iden- 
tified earlier when x* * 0. Thus the asymptotics pre- 

. . / • „ i.1 


The usual flame-sheet analysis determines the gra- 
dient on the unbumed side of the flame. 


Viuuu vaiuv* ^ W - i. L 

diets a turning point (quenching point) in the re- 
sponse. but only provides a first-order estimate of its 
location (T* = T*, S' 2 = rtf 2). and not a very ac- 
curate one at that. 

A refinement of this description can only come 
from a higher-order analysis, and in this connection, 
we rewrite equation 4 in the form 


Yf l Y v (x* + 0) = [2DT^- 2 e-^] 1/2 


= {S[pC p /aX] x/2 )^ = S 


where S is the adiabatic flame speed defined by the 
system. (The label () + is used to denote a quantity 
that is expressed in terms of dimensional variables.) 
Matching with the gradient of equation Al(b) then 
yields a formula for x*, 


j2Jne ' d /2 /erfc(x* /y2) = $ 


and from this, we can deduce variations of x* with 
$, where 


X* = (SpC p A- l x*)+ - x*5 


is a nondimensional flame-sheet location that is not 
scaled with a. x* decreases monotonically with in- 
creasing a, reaching the value 0 when 3 _2 (~a) = 
ti/ 2, a well-known result (e.g., Ref. [13]). For these 
solutions, the maximum temperature remains fixed 


at T a = 1 4- Tf. 


There is also a solution for which x* — 0 and the 
flame temperature T* is within 0(1/0) of T a . Then, 
to first order, 


x > 0 T = Tj + eric(x/j2) (A5) 
Within the reaction zone, the variables are f , £ where 
T = r a (l + 0~ l T a t) y X = 0~ l £ (A6) 


whence 


t ie = D0~ 2 e~ H/Ta te t 


t = t*. t, = 0 at i - 0 


Here t* defines the flame temperature, a quantity to 
be determined. 

A single integration leads to 


(f,) 2 = 2D0- 2 e-» n * [p(t*) - p(t)]. 


p(t) - e* (1 - t) 


r, f + (nt2)Der (T a - T)e^ r « 0. 


= erf(x/jz) 


Now the turning point can onlv be generated when 
at least the exponential tail of the temperature per- 
turbation within the reaction zone interacts with the 
symmetry plane x = 0. The natural assumption is 
that x* = O(ln0/0) when this occurs, corresponding 
to e’ k * = 0(1/0) for some k . Then we can approx- 
imate by 1, and that this is a poor approximation 
beyond the reaction zone is obviously of no concern. 
Equation A10 can then be integrated using the 
boundary conditions 


T s - 0, T = T* at s = 0 


to yield 


dT [h(-0/T+) - h(-6/T)Y 


= - Jn9D.s 


hit) * 1 0[E,(t) - t-'e 1 - t- 2 e‘] 


which imposes the requirement t < f* < 0, a range 
in which p is an increasing function. Matching be- 
tween the gradients defined by equations A5 and A8 
(x — > 0, t -* — o°) then leads to a formula for f*, 
namely. 


+ T a [E x (t) - t^e<) 

where E,(x) = / r _. dt r 'e* is the exponential inte- 
gral. Equation A12(a), with s = 1 (x = °°), T — Tf* 
is an implicit formula for T+{D)> and this function is 
also plotted in Fig. 1 (curve A). It is a decent ap- 
proximation to the exact solution. 

If the problem implicit in Equations A5-A7 is ex- 
amined to the next order, in an attempt to calculate 
a better approximation to the flame-temperature for 
the wall flame, that is, t* = f * \ + 0 _1 f*2. & ^ not 
difficult to show that t * i — » 00 a st^ l — > 0, with 
— » const., give or take logarithmic terms. (Second- 
order flame-sheet theory is seldom required in our 
subject, but a detailed example is presented in Ret. 
[14].) Thus the erstwhile 1/0 2 term in T* is corn- 
parable to the 1/0 term if = O(l/j0)> corre 
sponding to an 0( 1/0^0) deviation from T a . The 
turning point is characterized by this scale. 
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{ satisfy the boundary condition at s = 1, where 
;*)* a constant. 

now equate T as defined by A 16 with A18 
ag to all orders) so that A — — Tlj2Dq(t+) 

I V 2 ^ (t.) = (T„ - Tf) + Tie X 

k+ ri* |[i " - 111 

formula that describes variations of the flame tem- 
with the Damkohler number Within the 
-text of the autonomous problem, equation A19 
.Valid to all algebraic orders since an arbitrary num- 
ber of terms can be retained in q. 



^SjWhen only the first term in q(t) is retained [the 
i that is nominally 0( 1)], the error in f is 0( I/O), 
B compared to 1 jJO, and the turning point is cap- 
l Fig. 1 (curve Al). Retaining two terms in q(t) 
i a better approximation, also shown in Fig. 1 
re A2). 

-The term proportional to 6 ' 1 in A19 plays a crucial 
- in defining the turning point; for at all orders, 
m = 0 so that the integral is logarithmically sin- 
as f* — > 0. (That all terms in the expansion of 
F (equation A14) have a zero derivative at the origin 
^ns that the expansion of each q-dependent term 
■ A 19 for large 0 is uniformly valid in f*. further 


justification for the validity of replacing q by its first 
approximation.) Indeed, for small f* and only the 
first approximation for q , A19 can be approxi- 
mated by 

Tj j2Dq{t„) ~ (T. - T { ) + Tjtf' 1 X 

[ - Jz ln{ — f*) + 1.34992 . . .] (A20) 

see Fig. 1 (curve B). The derivative of A20 is ap- 
proximately 

(-1/2) jl fit* + (vji D)dD/dt+ 

~ -fo(0t+) (A21) 

which exphcitlv shows that in the immediate neigh- 
borhood of the turning point f* is 0(1/ JO). 

A final note. The formula A 16 may be used to 
provide a measure of the overall thickness of the 
reaction zone (bounded on the left by the symmetry 
plane). The displacement effect can be calculated by 
finding the value of f for which the linear variation 
of t , extrapolated to small £, defines a value oft equal 
to f*. Thus, 

- n(T a - 7>)- 1 x 

[ — J 2 . ln( — f») + 1.34992...] (A22) 

which is 0(ln 0 ) when t, is 0(1/ Jd). Recall the re- 
marks following equation A 10. 
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COMMENTS 


J. Chomiak, Chalmers University of Technology, Swe- 
den. Your edge-flames are often curved. I wonder what is 
the mechanism of the branching phenomenon and how the 
branching affects the speed of the flames? 

Authors Reply. The fresh/fresh counterflow edge-flame 
is horseshoe shaped, which is what one would expect. As 
for the edge -flame constructed from a deflagration with 
heat losses, it has only two choices: to be straight or curved, 
and in a two-dimensional asymmetric combustion field, it 
would be surprising if it were straight. The same can be 
said for our previous example of a premixed edge- flame, 
the one generated in a fresh/inert counterflow (Ref. [8] in 
the paper). The most obvious effect of edge curvature on 
the propagation speed occurs when the Lewis number (Le) 
is different from X, and we will report results for premixed 
flames in this case elsewhere. As we briefly note for dif- 
fusion edge-flames, small values of Le lead to enhanced 
reaction near the edge, and the edge-speed remains posi- 
tive ail the way to the i-D quenching point. Indeed, edge 
advancement is possible beyond the 1-D quenching point: 
the edge trails cellular structures, sublimit flames— obvi- 
ously, it can not trail a 1-D flame. This can also happen for 
the deflagration -with-heat-losses problem, even if Le = 1, 
for reasons that we discuss elsewhere. 


Robert Pitz, Vanderbilt University , USA. Do you expect 
substantial differences in your conclusions if more complex 
chemistry is employed in your analysis? 

Author’s Reply . No. The essential characteristics of un- 
bounded edge-flames — well defined waves that can travel 
with positive or negative speeds, depending on the Da- 
mktthler number — have their roots in the fact that the un- 
derlying 1-D problem has two stable solutions. Whether 
the 1-D configuration is a diffusion flame or a premixed 
flame is a mere detail, as is the chemistry. Premixed failure 
waves (edge-flames with negative speeds) have, we believe, 
been seen in the sublimit failure of upward propagating 
methane-air flames in a standard inflammability tube, as 
we note in (Ref. [8] in the paper). Ronney has recently 
been able to establish stationary premixed edge- flames in 
his laboratory, and his observations are consistent with the- 
ory. These are analog-computer simulations with full chem- 
istry. We can safely say that there are three, not two, basic 
flames in combustion: diffusion flames, deflagrations, and 
edge-flames. 
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